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Vertical probes with temperature sensing elements placed at fixed positions along their length 
are commonly used for measuring earth temperature variations with time. Mainly for structural reasons, 
the probes are comprised of materials whose thermal properties are not the same as those of the sur- 
rounding earth, so that the temperatures as measured at a given time by the probe are not the same 
as that for the undisturbed earth. A mathematical analysis for steady periodic, two-dimensional heat 
flow in a two body composite has been made to determine the probe conduction errors in ground tem- 
perature measurements. Several examples are given to show the relative magnitude of probe con- 
duction errors. 
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1 . Introduction 

The diurnal and annual variation of the tempera- 
tures in earth or soil at varying depths beneath the 
ground surface has a considerable significance in a 
wide variety of applications. These include the de- 
termination of thermal properties of undisturbed soils 
from temperature measurements, the depth at which 
water and sewer pipes have to be buried to prevent 
damage by freezing, the effect of temperature varia- 
tions on heat losses from buried utility pipes and the 
design of house foundations, and the need to know the 
temperature as a function of time and depth adjacent 
to unoccupied underground cavities in order to eval- 
uate heat transfer in the earth from time of occupancy. 

The accurate measurement of temperature varia- 
tions demands that some attention be given to the 
method of measurement. Some of the considerations 
necessary for the method are that the earth be rela- 
tively undisturbed by the insertion of temperature 
sensing elements (usually thermocouples, resistance 
thermometers or thermistors), and that the depth of 
the element below the surface of the earth be known 
quite accurately. These considerations are usually 
met by attaching thermocouples securely at fixed 
positions along the length of a rigid probe (pipe or 
tube) before inserting the probe into the earth. The 
positions or depths below the surface of the thermo- 
couples are then known quite accurately. Other con- 
siderations for the accurate measurement of 
temperature are the protection of the temperature 
sensing elements from corrosion and stray electric 
fields. 



The probe which is inserted into the undisturbed 
earth must be sufficiently strong and rigid that it not 
be deformed or broken upon insertion. This structural 
requirement may not be compatible with the re- 
quirement that the thermal properties of the probe 
material be similar to the thermal properties of the 
surrounding earth. Under the influence of temperature 
variations at the surface of the earth, appreciable 
differences in thermal properties will lead to errors 
in the measured variation of temperature. 

The purpose of this paper is to show the effect on 
temperature variations at depths below the surface 
by placing a vertical probe in the earth. 



2. Mathematical Analysis 



This analysis of steady, periodic, two-dimensional 
heat flow assumes a simple model such as shown in 
figure 1, where a sine-wave temperature variation is 
applied to the surface, 2 = 0, and at a sufficient depth, 
z=/, below the surface, the temperature is invariant 
with time. Region 1 comprises the earth with a thermal 
conductivity k\ and thermal diffusivity a.\. Region 2 
is comprised of the probe material with a thermal 
conductivity \ 2 and thermal diffusivity a 2 , and is in 
perfect thermal contact with region 1. Other more 
complex models and boundary conditions can be as- 
sumed using a similar analysis. 

For a cylindrical geometry with the axis of the sys- 
tem at the center of region 2 (fig. 1), the partial dif- 
ferential equation for heat conduction in two regions is 
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FIGURE 1. Model for steady periodic heat flow in earth {region 1) 
shown with probe {region 2) of radius, r=a. 

At surface 2=0 temperature varies sinusoidally and at z=l the temperature is invariant 
with time. 
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where y 2 n = q\+ (nn/l) 2 , ft = q\ + (nir\l)\ U and K 
are zero order modified Bessel functions of the first 
and second kind respectively, and A n and B n are con- 
stants to be determined from (5) and (6). Substitution 

of (7) and (8) at r= a in (5), yields 
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Multiplying both sides by sin kirz {k =1,2. . .) and 
integrating from z = to /, gives 
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where v is the temperature potential measured from 
the temperature at z = /, a is the thermal diffusivity 
and subscript i denotes either region 1 or 2. There are 
several methods by which this problem may be 
solved and the author has found it convenient to use 
the Laplace transform method. The Laplace trans- 
form of (1) is 



d 2 v j 1 dvi d 2 Vj 
dr 2 r dr dz 2 
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where q 2 = p/cxi and p is the parameter of the Laplace 
transform. 

For region 1, the radial temperature gradient is 
assumed to approach zero for very large values of 
the radius, r. The boundary conditions and con- 
tinuity conditions at the interface of the two regions 
(fig. 1) with their respective Laplace transforms are 
as follows: 
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where d= I 2 1 a 2 n 2 7r 2 and c= l 2 \a x n l TT 2 . Differentiation of 
(7) and (8) with respect to rand substitution in (6) yields 
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Substitution of (10) in (9) gives A n , which when sub- 
stituted in (7) gives at the interface, r=a, 
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where (x) = 2ttIT, and T is the period of a cycle. Solu- 
tions of (2) that satisfy (3) and (4) for the two regions are 



Vi=- 



o) sinh q\(l — z) 
p 2J r(j) 2 sinhgi/ 



AnKo(ynr) . jittz 



Ml+W) 1/2 /i(L)£oC/)J 



r(p)=(i+pc)(i+prf)|i + 



L = ^p(l+p^ 



To obtain vi from V\, the Inversion Theorem for the 
Laplace transformation is used, where for determining 
the steady periodic portion only the simple poles at 
p = ±ia) need be considered. The residues at the 
poles yield 
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where tff = vTrl{a.\Tll 2 ) 

X= sinh i// cos i// 

F= sin i// cosh i// 

/? = sinh i// (l — z/l) cos i/j (1-z/Z) 

5= sin (1-z/l) cosh i// (1-z/l) 

W(±iu))=P±iQ 

The real and imaginary parts of JF are determined 
using complex arithmetic where the real and imaginary 
parts of the Bessel functions of complex argument are 
given by 

In(pe l *)= Mnip, <p) + iN„[p, if) 

Kn(pe^)=G n {p, (pHiH n {p, if). 

A computer subroutine has been programmed to give 
values for M„, /V„, G n , and H„ for n= and 1. Numeri- 
cal values for the temperature potential at the interface 
between the probe and the earth have been calculated 
from (12) using a Univac 1108 digital computer. 

3. Numerical Solutions 

If the thermal diffusivities of the earth and the probe 
material are equal, the last term on the right side of 
(12) goes to zero and temperature variations with time 
and depth from the surface will be the same at corre- 
sponding depths anywhere in the earth. 

For numerical solutions, there are numerous pos- 
sible variations of the parameters. For this reason, only 
the effects of variations in the thermal properties of the 
probe material (region 2) will be studied while the other 
parameters will be held constant. For a diurnal tem- 
perature variation, let 

r=24 hr (86400 s) 

/ = 5 ft (1.524 m) 



A, =0.21 BTU hr ' ft ' (deg F)" 1 (0.3633 w m 1 K- 1 ) 

a, = 0.01 ft 2 hr 1 (2.58 X 10 < m 2 s ' ). 

Because the probe may not be comprised of a homo- 
geneous material and in most cases will be a filled 
tube, its thermal properties are difficult to specify. 
For this reason, four sets of "effective" thermal prop- 
erties were investigated characterizing undisturbed 
earth and low, medium and high thermal conductivity 
probe materials. The thermal properties for case 2 may 
approximate those for wood or bake lite, for case 3, 
stone or masonry and for case 4, stainless steel. 

Numerical solutions of V\ against time (12) for the 
four cases are shown in figure 2. Curves are shown for 
the 2 in (0.0508 m), 6 in (0.1524 m) and 12 in (0.3048 m) 
depths and for only a 12 hr portion of the 24 hr cycle. 
Since the amplitude of the temperature potential was 
assumed to be unity (3), the actual temperature 
0(<7, z, t) at the radius r=a as a function of time and 
position may be found from the ratio 
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where 6\ is the temperature at z = I and 02 is the ampli- 
tude of the temperature variation at z = 0. 

In case 1, the probe material has the same proper- 
ties as the surrounding earth and therefore the temper- 
ature variation is the same as that for undisturbed 
earth under the boundary conditions cited. Tempera- 
ture variations for cases 2 and 3 give negligibly small 
deviation from the undisturbed earth temperature at 
the 12 in depth. Case 4 representing a highly conduc- 
tive probe shows large deviations from the undisturbed 
earth temperature at these depths. With this probe 
there is a considerable phase angle deviation which 
increases with increase in depth below the surface. 

Figures 3a, 3b, and 3c show the temperature differ- 
ences between the temperature potential at r=a and 
that for the undisturbed earth for case 2, 3, and 4 as a 
function of time and position below the surface z=0. 
Letting A be the temperature difference as described 
above, then the actual temperature difference between 
the temperature of the probe at r=a and that for the 
undisturbed earth is 
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BTUhr^ft" 1 (degF) 1 


wm-'K" 1 


ft 2 !!!- 1 


m 2 s ! 


Case 1 0.21 
(Undisturbed earth) 
Case 2 0.08 
("Low" eonduetivity) 
Case 3 0.8 
("Medium" conductivity) 
Case 4 10. 
("High" conductivity) 


0.3633 
.1384 
1.384 

17.3 


0.01 
.0044 
.03 
.4 


2.58 X10 7 
1.135 X10 7 
7.742X10 7 
1.032 XIO 5 
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A„ = 2 A. 

Assuming a 10 °F amplitude for the diurnal variation, 
then for case 4 at t = 6 hr the actual temperature 
difference is 2.26 °F for the 2 in depth, 3.06 °F for the 
4 in depth, 2.66 °F for the 6 in depth and 0.45 °F for 
12 in depth. 



the actual surface temperature variation for a weather 
cycle is given by the trigonometric series 

K 

^=2 B m sin (a> m t — <p m ) 



4. Discussion 

The mathematical development given in a previous 
section considers a pure sine-wave temperature 
variation at the surface of the earth, whereas in nature 
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FIGURE 2. Temperature variation with time and depth below the FIGURE 3. Temperature difference between temperature at r= a 
surface at the radius r= a {fig. l)for the four cases, assuming unit and that for undisturbed earth for (a) case 2, (b) case 3 and (c) case 4 

amplitude in the temperature variation at the depth z = 0. as a function of time and position below the surface z = 0. 
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where k can be considered to be a finite number. This 
relationship can also be satisfied in the mathematical 
development by its substitution for (3). This was not 
considered necessary for the purpose of this paper, 
which is to present an analytical method and to show 
by numerical results the relative magnitude of probe 
conduction errors on a diurnal basis. If an annual cycle 
is to be considered, the depth /, at which the tempera- 
ture is assumed to be invariant with time should be 
much greater than the 5 ft used for the diurnal cycle. 
A study of the transient case was made by Brown. 1 
He performed tests which produced a monotonically 
increasing temperature change with time at the various 
depths below r the surface. For the usual periodic 



1 W. G. Brown. Thermal model tests for probe conduction errors in ground temperature 
measurement (Res. Paper 194, Division of Building Research, National Research Council, 
Canada, Sept. 1963). 



weather cycle (diurnal and annual), the present 
analysis permits estimation of the magnitude of the 
harmonically varying probe conduction errors repre- 
sented by the differences in temperature when employ- 
ing probes of different thermal properties (cases 2, 3, 
and 4) compared to that of undisturbed earth (case 1). 
The impetus for this paper is the determination of 
the thermal properties of undisturbed soils or earth 
from ground temperature measurements. This work is 
currently sponsored by the Office of Civil Defense, 
Department of the Army, and consists of this analysis, 
laboratory tests of probe configurations under time- 
varying temperature conditions, and actual field 
measurements now underway. 

(Paper 72C4-281) 
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